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The Lefschetz thimble method, i.e., the integration along the steepest descent cycles, 
is an idea to evade the sign problem by complexifying the theory. We discuss that such 
steepest descent cycles can be identified as ground-state wave-functions of a supersym¬ 
metric Hamilton dynamics, which is described with a framework akin to the complex 
Langevin method. We numerically construct the wave-functions on a grid using a toy 
model and confirm their well-localized behavior. 

Subject Index A22, A24, B15 


1. Introduction. The first-principle approach to solve the theory is an ultimate goal of 
theoretical investigations. The quantum Monte Carlo simulation for the functional integral 
has been a powerful ab initio technique to reveal non-perturbative features of various phys¬ 
ical systems. Successful applications include the lattice-QCD (Quantum Chromodynamics) 
simulation, the lattice Hubbard model, the path integral representation of spin systems, etc. 

The Monte-Carlo algorithm is based on importance sampling, so it is demanded that the 
integrand should be positive semi-definite, i.e., e'^ > 0 where S is the classical action. This 
positivity condition is often violated in systems of our interest and then we can no longer 
rely on Monte Carlo simulations [1, 2]. In QCD a finite baryon or quark density introduces 
a mixture of Hermitian and anti-Hermitian terms in the action, and then e'^ acquires a 
complex phase. In repulsive Hubbard model away from half-filling or generally in fermionic 
systems with spin imbalance, the sign of the integrand may fluctuate. It is also a notorious 
problem of the complex phase that appears from the Berry curvature in the path integral 
representation of spin systems and this phase cannot be removed for frustrated situations 
such as the XY model on the Kagome lattice. 

Moreover, to approach real-time quantum phenomena, e*^ is an oscillating function by 
definition and the sign problem is unavoidable. Although the Monte Carlo simulation is 
useful to compute physical observables in equilibrium in the imaginary-time formalism, the 
analytical continuation is necessary to access the real-time information. In general, however, 
the analytical continuation is a quite pricey procedure, and some additional information on 
the system such as the pole and the branch-cut structures would be necessary. 

Many ideas have been proposed so far to overcome the sign problem and, unfortunately, 
applicability of each method has been severely limited. Recently new techniques to complex¬ 
ify the theory are attracting more and more theoretical interest, which includes the path 
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integral on Lefschetz thimbles [3-14] and the complex Langevin approach [15-21]. Except 
for several formal arguments, theoretical foundations for the complex Langevin method are 
not fully established, and not much is known about its reliability [16-18]. In the context of 
real-time quantum systems, the numerical simulation works for some initial density matri¬ 
ces [19, 20]; however, it is recently reported in Ref. [21] that the real-time anharmonic 
oscillator at zero temperature converges to a wrong answer with unphysical width. 

In contrast, the Lefschetz-thimble method has a solid mathematical foundation at least 
for finitely multiple integrals [3-5]; however, its practical applicability is still in the devel¬ 
oping stage. This method decomposes the original integration cycle into several steepest 
descent ones, called Lefschetz thimbles, using complexified field variables. Picking up a sin¬ 
gle Lefschetz thimble, one can employ importance sampling [6-8] and one can evade the sign 
problem since the oscillatory factor in totally disappears on each Lefschetz thimble. On 
the other hand, zero-dimensional model studies have exemplified the importance of struc¬ 
tures of multiple Lefschetz thimbles [10-13]. For further applications it is needed to deepen 
our understanding on more aspects of the Lefschetz thimbles. 

The purpose of this Letter is to shed new light on the Lefschetz thimble method in a 
form of the Hamilton dynamics, which was first elucidated in Ref. [5]. In this reformulation, 
the Lefschetz thimbles can be identified as ground-state wave-functions of a supersymmetric 
topological quantum system. After reviewing this modified Lefschetz thimble method, for a 
quartic potential problem at zero dimension, we solve the Hamilton dynamics concretely to 
find the corresponding wave-functions. Based on our analytical and numerical observations, 
we discuss advantages of this reformulation for the numerical computation. 


2. Lefschetz thimble and SUSY quantum meehanies. Let us consider an A^-dimensional 
real integral as a “quantum field theory” defined by a (complex) classical action S{x). In 
this theory our goal is to compute an expectation value of an “observable” 0(x) defined by 

/ OO 

d^x e^(^) 0(x) , (1) 

-OO 


where x = (x^^\x^^\ ..., x^^^) G and the normalization Af is chosen such that (1) = 1. 
The starting point in our discussion is to reformulate this theory in an equivalent and more 
treatable way using a complexified representation: 


(O) = I d^zd^z P{z,z)0{z) 


( 2 ) 


Here -|- ix^*^ and — ix^*^ with x^*\x 2 ^ G M and / dzdz represents the 

integration over the whole complex plane; i.e. dxi ff^dx 2 - The choice of the gener¬ 
alized weight function P(z, z) may not be unique. Indeed, a trivial example is P(z, z) = 
j\fQS{z) — z^'^'l). At the cost of complexifying the variables, nevertheless, it is often 

the case that P{z^ z) could be endowed with more desirable properties for analytical and 
numerical computation than the original 

A clear criterion to simplify the integral is to find P{z, z) such that the phase oscillation can 
be as much suppressed along integration paths as possible, while in the complex Langevin 
method P{z, z) is optimized to become a real probability. To suppress the phase oscillation, 
let us pick up a saddle point Za- satisfying S'{za) = 0. The steepest descent cycle or the 
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Lefschetz thimble of the saddle point is defined with a fictitious time t as 


Jo 



xi(0) + ix2(0) 


dxj*^ {t) 
dt 


dReS 


lim {xi{t) + ix 2 {t)) 

t^ — OO 



( 3 ) 


This is a multi-dimensional generalization of the steepest descent path in complex analysis, 
which we will refer to as the downward path. The original integration path on the real axis 
in Eq. (1) can be deformed as a sum of contributions on weighted with an integer rria', 
i.e., 4 m d^x = nia- fj z. In mathematics it is established how to determine from 
the intersection pattern between the steepest ascent {upward) path from z^j and the original 
integration path [3-5]. It is important to note that ImS is a constant on each Lefschetz 
thimble for the application to the sign problem [6, 8]. 

In the following, let us restrict ourselves to = 1 for simplicity, because the generalization 
is straightforward. So far, the Lefschetz thimble is constructed as a line, and let us find a 
two-dimensional smooth distribution P{z,z) according to Ref. [5]. For that purpose, we 
define the “delta-functional one-form” 5{J(j) supported on the Lefschetz thimble so that 



5{J^) A 0(z)e^(^4z. 


(4) 


For instance, (I(M) = 6{y)dy. Such delta-functional forms 5{Ja) (on a Kahler manifold) have 
a path-integral expression from the supersymmetric quantum mechanics [22-24] (see also 
Secs. 2.8 and 4 of Ref. [5] for more details in this context). Integration (4) can be represented 
as 


{O) = M J 'D[x,p,Tr,'ip] exp 


,' dxi (9Re S \ 


X exp 


/d , d^ReS , , 

dtTTii —5ij + 4 


dt 


dxidx^ 


O(z(0))e^("(°))(V'i+iV’2)(0) . (5) 


Here x,p are bosonic fields and tt, are fermonic ghost fields, and z{t) —)■ Zo- as t —)■ —oo. 
We should note that an integration in terms of z is promoted to the path integral on z{t) for 
t < 0, while the observable and the weight O{z{0)) exp S'(2:(0)) are functions of z{0) only. Let 
us outline how these two expressions (4) and (5) are equivalent [5, 22-24]. We first integrate 
out p{t) to get the Dirac delta function. 


Vp exp 



dReS \ 

dxi J 


/ dxi (9Re5\ 
\ dt dxi ) 


( 6 ) 


This delta function constrains the path integral on x{t) to a gradient-flow line defining 
Lefschetz thimbles. Since 2 ;(—oo) —)• z^, this path integral for t < 0 gives a delta-functional 
support on J^. However, the delta function produces an unwanted determinant factor. As is 
well-known, the path integral on ghost fields 7r(t), ^(t) for t < 0 can eliminate that factor as 


'DTr'D'ifj exp 


, , d , d'^ReS 

dtTTi \ —dij + 

dt oxiOXj 


4 


^ . d ^ d^ReS 

et( d^ 4 '+ 


(7) 


Now, we obtain an integration over surface variables x{0), V’(O)) and denote them by x, if). 
Locally, the Lefschetz thimble dla can be expressed as zeros of a certain function /, then we 
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can find that the path integral (5) eventually gives 


d^xd^Tp5{f) A 0{z)e^^'''> (V'i+iV’ 2 ) 

OXi 


5{f{x))df{x) A 0{z) dz, 


( 8 ) 


which is nothing but the local expression of the original integration (4). Going back to (5), 
this shows that the so-called residual sign problem comes from the fermionic surface term 
V'i(O) -|- iV' 2 ( 0 ) because one can identify 'ipi{0) = dxi as above. 

Importantly, with these added fields, pi,TTi,^i, the action is BRST exact under a transfor¬ 
mation; Sxi = dipi = 0, dvr* = —ipi, Spi = 0. By definition the nilpotency 5^ = 0 is obvious. 
Thanks to the boundary fermionic operator in (5), the surface term is BRST-closed so long 
as the observables are holomorphic. This makes a sharp contrast to the complex Langevin 
method that could also acquire the BRST symmetry but it is violated by the surface term. 
Because of the BRST symmetry we can add any BRST exact terms without changing the 
original integral, and it is useful to insert ^ f dtpf. In summary, the effective Lagrangian 
that describes the fictitious time evolution is given by the following topological theory: 


-heff — 


Ei 2 , ■ f dx 

-2ft 


j (9Re S \ 

“I j 


= -5 


.Ei , dxi 
TTi I l-Pi + ^ + 


dxi 
dReS 


dt 






d'^ReS 

dxidxn 




dt 


dxj 


(9) 


which is nothing but a Legendre transform of an effective Hamiltonian: 


Hes — ^ 

i 



1 

2 


fdReS^ ^dReS\ 

^ - Pi+Pi^ - 

\ OXi J 


E 


Id^ReSr.. ; 
2 


( 10 ) 


with [xi,pj] = iSij and {'n'i,'ipj } — Sij- The fermion number F = yfi'^i -|- TT 2 ip 2 is a conserved 
quantity of this Hamiltonian. After the time evolution from t = —00 only the ground state 
with the lowest energy eigenvalue remains, so that the generalized weight is given by 
P{z, z)dzdz = d>(z,z) Ae'^i^idz, where df(z,z) is the ground state wave-function and con¬ 
verges to in the limit Ei —)■ -|-0. Note that the weight factor exp 5 ( 2 ;) is necessary 

in this formula, since the wave function designates only the integration cycle Ja- We can 
further simplify this Hamilton problem by choosing e = Ei = £ 2 - Performing the conjugate 
transformation T = derivative terms are eliminated as 

-V- 

^ 2 

This describes supersymmetric quantum mechanics with the superpotential Re S [5]. 

Before applying this method to an interacting model, let us convince ourselves of 
perturbative correctness. For that purpose, we consider the simple Gaussian case: 

So(x) = -^x^, ( 12 ) 

where x G M is a one-component variable and uj G C. One can regard this Gaussian integral 
as an elementary building-block of perturbative quantum field theory; in a non-interacting 
theory of the scalar field cj), the action is decomposed into — ^[T — i(/c^ — w?)](j){—k)cl){h.) for 
each Fourier mode. In this case we can immediately hnd that the bosonic wave-function 
should be the ground state of harmonic oscillator with the ground state energy |a;|. The 


d^ReS 

dxidxj 


_TTi,'4)j 


( 11 ) 


^:fi = E 


nPi + WZ 


1 fdReS\ 


2i 


2e \ dx. 
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fermionic ground state energy — |w| cancels the bosonic energy thanks to supersymmetry. As 
a result, unoccupied fermions point the direction of the Lefschetz thimble; i.e., the super- 
symmetric vacuum belongs to the F = 1 sector. Let us see this in the simplest example, 
w G M and a; > 0. The fermionic part of the Hamiltonian is ^([7ri,'0i] — [^2,V’2]), and thus 
the 1 and 2 fermions are unoccupied and occupied, respectively. This leads to the fermionic 
ground state energy —a;, which cancels the bosonic ground state energy cj. Therefore, the 
unoccupied fermion is tangent to the Lefschetz thimble 17 = ]^. 

The final result of Po{z, z) for cj G C together with e'^ • is 


Po{z, z) = AAexp 


1^1 - , 1 / 2 I - - 2 \ ^2 

- — ZZ+—{0JZ +ojz)-—z 

2e 4e 2 


(13) 


which reproduces the original integral (1) with the action (12). To see this for a polynomial 
0(x) it is sufficient to require (z^) = Ijoj, which is nothing but the free propagator and can 
be explicitly confirmed with Eq. (13). In order for the exponentially fast convergence of Pqj 
the parameter e needs to be 0 < e < 2. We here emphasize that the theory is equivalent to 
the original (1) for 0 < Ve < 2 and the conventional Lefschetz thimble method is retrieved 
in the e —)■ 0 limit. Actually, in this limit of e —)■ 0, only a path of 2\oj\zz — loz‘^ — zoz^ (= 
{2lm{^/ujz)}‘^) = 0 contributes, which is nothing but a condition to guarantee Im5o(2:) = 
0 on the Lefschetz thimble. At finite e, this restriction is smeared and Pq may have a 
distribution around J with a width of order of e where a complex phase arises in general. 
The non-positivity of Pq is a big difference of this Lefschetz-thimble approach from the 
complex Langevin method, in which the distribution must be semi-positive definite. 

In the model with quartic interaction, S = —'^z^ — ^z‘^, the vacuum structure is drastically 
different from that of the Gaussian case (13). With A 7^ 0, there are three classical saddle 
points and the Morse index for each saddle point is 1. The Witten index is tr(—1)^ = —3, 
and thus there are three supersymmetric vacua in this interacting model for any e [25, 26]. 
In the path integral expression (5), we can distinguish these three vacua by specifying the 
boundary condition at t = —00, as we will discuss below in detail. 


3. Lefschetz thimbles at e —)• 0. We define our zero-dimensional model with the quartic 
interaction by 

S{x) = (14) 

and hereafter, we will specifically choose the model parameters as 

uj = l-i = \/2e“‘^/^, A = 1.5i . (15) 

This choice has been motivated by the application to real-time problems. As we already 
mentioned, Rew corresponds to the width T or e in the ie prescription, and Imw corre¬ 
sponds to —(/c^ — m^). Therefore, the parameters (15) represent a situation at k'^ = mf + T 
in quantum field theory. 

When e is small enough, the Hamilton dynamics should become equivalent to the con¬ 
ventional formulation of the integral on the Lefschetz thimble, which means that P{z,z) 
should have a peak along the Lefschetz thimble only. This makes the analytic treatment 
much accessible since we do not have to solve the quantum mechanical problem in this 
zero-dimensional toy model. To identify the Lefschetz thimble, we should integrate the flow 
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Re 2 Re 2 


Fig. 1 Changes of the Lefschetz thimbles for a; = 1 — 0.9i (left) and w = 1 — l.li (right) 
with A = 1.5i fixed. (Left) One of three thimbles as shown by the solid line contributes to the 
integral and two are to be dropped when w = 1 — 0.9i. (Right) All three thimbles contribute 
to the integral when a; = 1 — l.li. 


equation in Eq. (3) and find the upward and the downward paths. We show the numerical 
results in Fig. 1 for parameters slightly changed from Eq. (15). 

Eor our theory (14) three saddle points are located at zq = 0, z± = w/A. The latter 
is, for our choice of parameters (15), z± = ±(0.897 ±0.372i), and S(z±) = —1/3. We can 
see that a; = 1 — i is a critical value at which the destinations of the downward flows from 
the saddle points change drastically, that is known as the Stokes phenomenon, as is clear 
from two panels for cj = 1 — 0.9i (left) and a; = 1 — l.li (right) in Eig. 1. In general we can 
show that the Stokes phenomenon occurs at Im(a;^/A) = 0, and when A is pure imaginary, 
this condition gives arg (a;) = — vr/d. Importantly, not only the downward flows but also the 
upward flows change and the intersection number of the original and the upward paths 
changes accordingly [4, 5]. In the case with uj = 1 — 0.9i only one upward path from zq 
crosses the real axis as seen in the left of Fig. 1, and so the Lefschetz thimble going through 
Zq contributes to the integral. In the case with w = 1 — l.li, on the other hand, three upward 
paths from zq and z± all cross the real axis, and all three Lefschetz thimbles contribute to 
the integral. 

Keeping the potential application to the real-time physics in mind, we need deeper under¬ 
standing on the Stokes phenomenon. In fact, it occurs at = m? ± L and so the thimble 
structure may fluctuate depending on the frequency k^, while in Euclidean theory it never 
happens because k'^ — m? is always negative (except for an unstable potential with m? < 0 ). 
It should be an interesting future problem to clarify the treatment of the Stokes phenomenon 
in the real-time systems [ 10 ]. 

For the same model with a different set of parameters, the Stokes phenomenon has been 
discussed in Refs. [27, 28] and the probability distribution P{z,z) in the complex Langevin 
method has been numerically computed. The important insight obtained there is that P{z, z) 
in the complex Langevin method looks localized but has a power decay at large j^;], which 
causes a convergence problem. Hence, it would be an intriguing question how P{z, z) in the 
modified Lefschetz thimble method for e ^ 0 should behave especially at large \z\. 
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Wave-functions at finite e. As an application of the modified Lefschetz-thimble 
method with a regulator e, let us compute the wave-function for the e = ei = £2 case, 
from which P{z,z) is to be constructed immediately. We can readily find the eigenstate in 
terms of ifi for the last term in the Hamiltonian (11). By restricting ourselves to the F = 1 
sector, we can define the effective potential in a form of 2 x 2 matrix-valued function that 
amounts to 

/ 9^Re5/(9xf S/dxidx2 \ 

\ 9^Re5/clxi(9x2 —d‘^P.eS/dx\ ) ’ 




dReS y f dReS y 

dxi J I dx2 J 


and then the Hamiltonian is 


d' 


92 




2 


5x2 


(17) 


Let us solve the ground state of the above Lf'g at finite e, and we specifically adopt e = 1 
unless stated explicitly. 

When we solve the Hamilton dynamics, we should set initial conditions to select proper 
Lefschetz thimbles out. We choose the semiclassical ground state in the limit £ —)• -|-0 as our 
initial condition. Since 'I'' = and —)■ 5{J'a) in £ —)■ -|-0, ~ for 

small £. Let us consider the Lefschetz thimble around zq for instance, then its tangential 
direction at zq is given by X 2 = tan (7r/8)xi as seen also from Fig. 1. The initial wave-function 
is thus proportional to 


gRe5'(2;)/£^|-_ sin(7r/8)xi -I- cos(7r/8)x2) • (— sin(7r/8)dxi -|- cos(7r/8)dx2) . (18) 


The bosonic part is well localized at zq, which justifies the following choice of the initial 
wave function, = — 00 ) = 6{z — zq){— sin( 7 r/ 8 )dxi -|- cos( 7 r/ 8 )dx 2 ). Similarly, we can 

fix the initial wave function for z± as = 5{z — z±)(cos( 7 r/ 8 )dxi -|- sin( 7 r/ 8 )dx 2 ). 

For the numerical procedure we smeared the delta function in the initial wave function 
by a Gaussian as 5{z — Za) —)■ exp{—20[(xi — xicr)^ -|- (x 2 — X 2 a)‘^]}- Then we discretized xi 
and X 2 from —2.5 to -|-2.5 with dx = 5 x lO”^, We then numerically integrate = 

—using the Euler method with dt = 10“^ until the wave-function converges. The 
convergence is fast and stable and the wave-function hardly changes after t = 1 ~ 2. We also 
mention that we utilized the Crank-Nicolson algorithm to improve the numerical stability 
when we compute the Laplacian. 

With this prescription, we find three independent ground state wave-functions = 

‘'Mxi -|- ^ 2 ^ ‘'mX 2 , which is shown in Fig. 2. The figures (a,b) in Fig. 2 show two com¬ 
ponents of and the figures (c,d) show those of Since our system is symmetric 

under the reflection 2 ; 1 —)• — z, we can easily read from the figures of Remarkably, 

these ground-state wave functions are linearly independent from one another, which means 
that supersymmetry is unbroken. We have also verified that the Dyson-Schwinger equation, 

A(2;"‘)-I-a;(2:2) = (19) 

is satisfied within 1% accuracy. This clearly shows that the supersymmetric quantum 
mechanics provides a suitable framework to compute Lefschetz thimbles. 

Let us comment on the convergence of the wave-functions in the present modified Lefschetz 
thimble method. In the effective potential (16), the first term gives the dominant binding 
potential in our toy model for large \z\ because the first term is a polynomial up to the sixth 
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c 

u 





(c) 'I'j*'^+\xi,X2) 


Fig. 2 SUSY ground state corresponding to the saddle points 2^0 = 0 and z+ = 

^—ojjX with e = 1, w = 1 — i, and A = 10“^ + 1.5i. Denoting + 'I' 2 ^"^dx 2 , 

(a) and (b) represent while (c) and (d) represent 'I/'U+l. 





Fig. 3 SUSY ground state ^'(^ 0 ), which is obtained starting with the saddle point zo but 
with the initial relative weight equal to both components. 




order, while the second term is up to the quadratic order. Therefore, the convergence in the 
present numerical approach is quite improved and there is no problem of the power-decay 
unlike the complex Langevin method. One might think that the remaining is still 

oscillating, but this is not a problem practically. We have numerical verified that the effect 
of is very small in the region where the profile of is localized. 

We also checked the robustness of the numerical results against small variations of the 
center position and the smearing width in initial conditions. If the smearing width of the 
initial wave function is changed, the overall normalization would be changed naturally, but 
once we normalize the wave-functions in the same manner (in Fig. 2 we normalized them 
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<3 
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(b) ImP at a; = 1 — i, e = 0.2 


I 



Fig. 4 Weight functions P{xi,X 2 ) corresponding to the saddle point zq with e = 0.2. 
(a) and (b) show the result for a; = 1 — i, and A = 10“^ + 1.5i. (c) and (d) show the result 
for to = 1 + i, and A = 10“^ + 1.5i. 


as f dtci <lx 2 [{^'i^'"'^)‘^ + = 1), then we eventually get the same result. Such insen¬ 

sitivity implies that each wave-function is well localized and the overlap at the saddle point 
is small. However, the overlap at the saddle point is not completely zero. To see this effect, 
let us skew the relative weight of the initial condition so that the initial wave-functions are 
not orthogonal. In Fig. 3, we show and starting with the same relative weights, 

namely, 00 ) = S(z — zo)(dxi -|-dx 2 ) at the initial time, for the demonstration pur¬ 
pose. The result in Fig. 3 is given a natural interpretation as a superposition of three 

ground states shown in Fig. 2. Because of this overlap among wave functions, our prescrip¬ 
tion for the initial wave function needs further refinement in order to extract one Lefschetz 
thimble at e = 1. 

In order to see the importance of such refinements from another viewpoint, we checked 
behaviors of distributions P for small e. Let us set e = 0.2 with w = 1 — i and A = 1.5i, then 
we obtain Figs. 4 (a) and (b). Here, the normalization is given by f d^x P = 1. Although it 
is localized around the saddle point zq = 0, its shape is still far from a one-dimensional line 
shown in Fig. 1. In order for comparison, we show the result also for e = 0.2 with a; = 1 -|- i 
and A = 1.5i in Figs. 4 (c) and (d). For this case, P is well localized to a one-dimensional line, 
which is nothing but the Lefschetz thimble at this parameter. We numerically observed 
that the weight function P is well localized around a Lefschetz thimble if Stokes phenomena 
do not happen and e is sufficiently small. Expectation values of operators, such as ( 2 ^), 
also give correct numbers within the numerical accuracy with such parameters. However, 
as e becomes larger, wave functions spread as two-dimensional distributions, and expecta¬ 
tion values of operators do not necessarily give correct values. Since the Dyson-Schwinger 
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equation (19) is satisfied, this problem must come from superpositions of the wave function 
with other ground states. 

It is an important future study to clarify dependence on this initial condition more system¬ 
atically, in order to take an appropriate linear combination of these wave-functions giving 
each Lefschetz thimble. This would be a key step to study how the Stokes phenomenon is real¬ 
ized at £ = 1. It also opens a new possibility to take into account multiple Lefschetz thimbles, 
since ground-state wave functions can be superposed by changing the initial condition. 


5. Discussions and Conclusions. In this Letter, the supersymmetric reformulation of the 
Lefschetz thimble integration was studied and its practical computation was discussed. The 
Lefschetz thimbles are now regarded as the ground sates of a supersymmetric Hamiltonian. 
Our computational scheme for the Hamiltonian system is essentially the same with that for 
the Fokker-Planck equation in the complex Langevin method. Those supersymmetric wave- 
functions are numerically computed for a zero-dimensional toy model with the classical 
action 5 = -|- Sint where 5o = and S'int = —^x‘^. Since the Morse indices of all the 

saddle points are the same, the number of saddle points must be the same as that of linearly 
independent ground states. 

For the Gaussian model with only, the above-mentioned situation is clearly true, which 
is explicitly checked by computing the wave-function analytically and constructing one super- 
symmetric ground sate. From this almost trivial example, one can learn an important lesson 
about the two-dimensional smooth distribution Po(z,z). In the “semi-classical” limit of 
e —7- 0, we recover a delta functional support along the Lefschetz thimble in the original 
formulation. For non-zero e, the phase oscillation arises away from the Lefschetz thimble, 
and the formulation nevertheless reproduces the correct expectation value. 

For the interacting model with 5int, we performed numerical computations for the super- 
symmetric quantum mechanics, and confirmed the existence of three linearly independent 
ground states by restricting ourselves to the F = 1 sector. Since the wave-function in the 
F = 1 sector consists of two components, we must set them in the proper initial condi¬ 
tions. We started from a localized wave-function in the vicinity of each saddle point, and 
our numerical computation shows a remarkable stability under small modifications on the 
initial conditions. This reflects the fact that all the saddle points are attractive unlike the 
complex Langevin method. At the same time we also found substantial dependence on the 
initial relative weight of these two components. This clearly indicates that we need a careful 
refinement of the initial condition to use the modified Lefschetz thimble method for numer¬ 
ical simulations, but it also opens a new possibility for some convenient scheme to take into 
account multiple Lefschetz thimbles. 

In the case of the Fokker-Planck system, the ground state is often uniquely determined and 
the initial condition dependence does not appear as in the case in ordinary quantum mechan¬ 
ics without any special symmetry. The Fokker-Planck operator of the complex Langevin 
method can be written as a functional integral in a similar manner, and it also shows the 
same type of the BRST symmetry. There are, however, several important differences in these 
two formalisms: First of all, the diffusion term in the complex Langevin equation does not 
give the gradient flow because the sign is different. Because of this difference, the BRST 
invariance cannot be promoted to the supersymmetric system satisfying 2ff = {Q,Q} with 
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supercharges Q, Q. Moreover, the fermionic sector should be restricted to F = 2 instead of 
F = 1. Therefore, it is not straightforward to relate these two formalisms yet. If one could 
establish a firm relation between two methods, it would be a great progress and the present 
modified Lefschetz thimble could provide us with an opportunity for a fully complementary 
approach to the sign problem together with the complex Langevin method. 
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